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We demonstrate the effect of bandstructure on the superconducting properties of Pb by calcu- 
lating the strong-coupling features in the optical conductivity, cr(aj), due to the electron-phonon 
interaction. The importance of momentum dependence in the calculation of the properties of su- 
perconductors has previously been raised for MgB2^'^. Pb resembles MgB2 in that it is a two band 
■ superconductor in which the bands' contributions to the Fermi surface have very different topologies. 

' We calculate a-{uj) by calculating a memory function^ which has been recently used to analyze cr(tj) 

of Bi2Sr2CaCu208+/. In our calculations the two components of the Fermi surface are described 
by parameterizations of de Haas-van Alphen data. We use a phonon spectrum which is a fit to 
C ' neutron scattering data. By including the momentum dependence of the Fermi surface good agree- 

ment is found with the experimentally determined strong-coupling features which can be described 
by a broad peak at around 4.5 meV and a narrower higher peak around 8 meV of equal height. 
5-H ' The calculated features are found to be dominated by scattering between states within the third 

band. By contrast scattering between states in the second band leads to strong-coupling features 
^ , in which the height of the high energy peak is reduced by ~ 50% compared to that of the low 

energy peak. This result is similar to that in the conventional isotropic (momentum independent) 
treatment of superconductivity. Our results show that it is important to use realistic models of the 
bandstructure and phonons, and to avoid using momentum averaged quantities, in calculations in 
order to get quantitatively accurate results. LA-UR-06-4303 
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PACS numbers: 74.25.Nf, 74.25. Gz, 74.25. Jb 



I. INTRODUCTION 



^ , Suhl et al.^ originally introduced a two-band model for superconductors as an extension of the original BCS 
'sj" ' calculation in which each band has its own superconducting gap. This two-band model was proposed to explain fine 
OO , structure in the infrared absorption of superconductors with s and d electrons. Many years later Allen® proposed a 
formalism to extend the Eliashberg treatment of superconductivity to include multiband Fermi surfaces and Fermi 
surfaces with complicated geometry by expanding order parameters, pairing interactions and scattering mechanism in 
terms of Fermi surface harmonics. This formalism was employed by Golubov and Mazin '' to investigate the effect of 
' impurity scattering on the value of the critical temperature in two-band model superconductors. Kresin and Wolf^ also 
" ^ . discussed two band superconductivity in the context of the high Tc materials in which superconductivity in one band is 
' induced in a second band by interband scattering. These investigations were not extended to get quantitative results 
for specific materials. Quantitative investigations of superconductivity in real multiband systems have only been 
I ■ developed in recent years with the availability of the necessary computational resources. Renewed interest in this area 
■ O [ is driven by experiment. In 2001 Nagamatsu et al.^ discovered that MgB2 becomes a superconductor at temperatures 
^ ■ below Tc — 39 K. This value of Tc is comparable to that of Lai_a;Sr2.Cu04, the first of the cuprate superconductors 
to be discovered. This discovery initiated a number of theoretical and experimental studies on the nature of the 
superconducting state in this material. MgB2 shows a significant isotope effect, a tell-tale sign of the involvement of 
phonons in the process^". Several tunneling experiments found very different values for the superconducting gap^^'^^, 
which led Liu et al. to suggest that the superconducting state of MgB2 exhibits multiple gaps. First principle 
5^ ' calculations"'^^ found essentially two distinct conduction bands, and it is this two-band nature of the MgB2 system 
I that is responsible for two different gaps. The calculated Fermi surface was subsequently found to be in agreement with 
angle resolved photoemission (ARPES)^^ and de Haas-van Alphen^® measurements. Choi et al. ^'^ calculated the gap 
function and Tc using the anisotropic Eliashberg formalism^^ which incorporates the momentum dependence of the 
phonons, electronic bandstructure and the electron-phonon interaction. Choi et al. used the ab initio pseudopotential 
density-functional method to calculate these. They pointed to the strong variation in the electron-phonon interaction 
on and between the Fermi surfaces of the two bands in which the interband scattering is much weaker than intraband 
scattering because of the different symmetry of bands in MgB2."''^ The strong dependence of the electron-phonon 
coupling on Fermi surface states in MgB2 leads to a wide range of values for the superconducting gap. Choi et 
al. emphasized the importance of including the momentum dependence of these quantities rather than using the 
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isotropic version of the Eliashberg equations for quantitatively accurate results. Here we show that there is also a 
strong dependence on bandstructure in the optical conductivity of superconducting lead. Pb is effectively a two band 
material in which the electron-phonon interband scattering is comparable to the intra-band scattering so that the 
superconducting gaps on the two bands are almost equal, in contrast to MgB2. Tc ^ 7.9 K, is also much lower in 
Pb than in MgB2. In order to investigate the importance of including the different character of the two bands which 
cross the Fermi level in Pb we calculate the strong-coupling features which appear in optical conductivity. 

Optical conductivity can be thought of as having two contributions. The first one is associated with the anomalous 
skin effect, whereas the second is due to "bulk" properties of the metal. The anomalous skin effect depends on the 
nature of the scattering of electrons from the surface and arises from the breakdown of momentum conservation at the 
surface. This contribution was initially investigated theoretically by a number of authors in the normal state^^'"'^^'^'^. 
In superconductors it was studied by Pippard using the phenomenological London model^" and by Mattis and Bardeen 
on the basis of the BCS theory^^. The second contribution, the "bulk" part, comes from scattering of the incident 
photons on phonons, impurities, or other forms of scattering centers, within a distance comparable to the skin depth 
from the surface. We will focus on the Holstein process, the simultaneous absorption and emission of phonons by an 
electron, and ignore impurity and electron-electron scattering in our calculations using a memory function approach. 
This has previously been investigated by AUen^^. 

Farnworth and Timusk^'^'^^ studied the strong-coupling effects associated with the electron-phonon interaction 
by infrared absorption. They measured the frequency dependence of the light emitted from a cavity whose walls 
contain Pb. They analyzed their data by taking the derivative of the ratio of the signal from the cavity in the normal 
and the superconducting states which leads to a frequency dependence given by the difference of the derivatives 
of the absorptivity in the superconducting and normal states, As{Ld) and A]s[{uj). This procedure has the effect 
of subtracting out the anomalous skin effect contribution and of enhancing the difference between the normal and 
the superconducting states. This results in a transport version of the a^F function which arises in the Eliashberg 
formalism. 

Farnworth and Timusk's data, shown in Fig. 1, exhibits features associated with the phonon density of states at B, 
C, D, E, F, and G. In addition the data has a feature, A, at w = 4 Aq, where Aq is the value of the superconducting gap. 
This feature could be thought of arising either from the phonon-mediated interaction between the superconducting 
quasi-particles or as an effect due to the change in the phonon self-energy in the superconducting state. The leading 
contribution to this feature is then fourth order in the strength of the electron-phonon coupling. The data of Farnworth 
and Timusk on thin films reveals a splitting of this feature into two slightly different but distinct energies indicating 
that we are dealing with two conduction bands in Pb^"^. 

We use the memory function formalism which was developed by Gotze and Wolfic'^ to calculate the contributions to 
<j{uj) due to scattering events from phonons, impurities, and spin fluctuations. The assumption of this approach is that 
the memory function, M{uj), can be expressed as a power series in the electron-phonon interaction strength, which 
makes this approach in principle a weak coupling approach. The advantage of this approach, compared to summing 
ladder diagrams, is that it is much easier to include the details of the Fermi surface which leads to improved quanti- 
tative agreement with experiment for Pb. Previously Tomlinson and Carbotte^^ have investigated the temperature 
dependence of the resistivity using the same experimental data as we use for Pb with different approximate solutions 
to the Boltzmann equation. Leung et al. have also taken a similar approach for the resistivity of aluminum. 

Recently Hwang et al."* determined M{uj), referring to it as the "optical single-particle self-energy" , in their analysis 
of <j{ijj) in Bi2Sr2CaCu208+5 and compared it with a single-particle self-energy derived from ARPES data. The 
interpretation of (t{lo) data in the high Tc materials, as well as other types of data, is considerably more difficult than 
in conventional electron-phonon superconductors. In high Tc superconductors strong correlation effects are evident 
in phenomena, such as the pseudogap, and are thought to be responsible for the high value of Tc- These correlation 
effects and the absence of a consensus on a suitable microscopic model make inversion of data and analysis based on 
weak coupling treatments problematic. Schachinger et al.^^ have recently investigated the inversion of (j{io) data for 
conventional and high Tc superconductors using singular value decomposition and maximum entropy techniques. 

In our calculations we use fits to the experimentally determined Fermi surface determined from de Haas-van Alphcn 
data and phonons from neutron scattering data, rather than calculate these as Choi et al. did in the case of MgB2 . 
In Pb there are two 6s and two 6p electrons per atom leading to four bands. De Haas-van Alphen measurements 
and a subsequent analysis of these by Anderson and Gold^^ show that in Pb the lowest lying band is completely 
inside the Fermi surface, and that the second and third bands both cut across the Fermi energy, whereas the fourth 
band lies above the Fermi energy. The component of the Fermi surface due to the second band has the topology of a 
sphere while that due to the third band is multiply connected and has the form of cylinders along each edge of the 
first Brillouin zone. The two band nature of Pb is due to the second and third bands. The phonons are described by 
Cowley's fit^^ to neutron scattering data. 

In sections II and HI we discuss our treatment of the charge carriers and the electron-phonon coupling. In section 
IV we derive an expression from the electron-phonon contribution to the memory function in the superconducting 
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FIG. 1: Reproduced from reference 24 with permission: Comparison of a F numerically inverted from infrared absorption 
(solid line) and tunneling (dashed line) 



state and an expression for the absorptivity in terms of the memory function. After calculating the derivative of 
the difference of the calculated absorptivities in the normal and superconducting states we compare the result with 
the data of Farnworth and Timusk^^'^^. The comparison demonstrates the importance of scattering involving the 
third band and resolves the discrepancy between the experimentally determined a^F and that calculated ignoring 
momentum conservation and approximating the Fermi surface with a spherc^^. 



II. TWO-BAND CHARACTER 



Anderson and Gold^^ carried out de Haas- van Alphen measurements on Pb and their subsequent analysis is based 
on the orthogonalized plane wave method which is fit to dc Haas-van Alphen data of the Fermi surface. The Fermi 
energy, ep, the Fourier components of the pscudopotential parameters, Vinand V200J f^nd a spin-orbit coupling 
parameter A are regarded as fitting parameters. We used the following values: ep ~ 9.765 eV, Vm = —1.142 eV, 
^200 = —0.530 eV, and A = 1.306 eV. Anderson and Gold's parameterization yields the dispersion relations for the 
4 electron bands. The Fermi surface cuts through the second and third bands, which are shown in Figs. 2 and 3. 
Comparing these figures with the form of the Fermi surface calculated for an empty face center cubic lattice'^^, one 
sees that the effect of interactions among the electrons is to remove the sharp edges in the second band and increase 
the density of electrons in the third band. Electron interactions also now ensure that the fourth band is empty at low 
temperatures whereas there is a small contribution to the electron density from the fourth band in the empty lattice 
calculation. By using this experimentally determined Fermi surface we have included the effects of electron-electron 
scattering in the dispersion of the second and third band carriers. 

In order to capture the strong-coupling features, the contributions from scattering among the different components of 
the Fermi surface have to be evaluated very accurately. These contributions are given by six-dimensional integral and 
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require a large number of mesh points. The Anderson and Gold model requires the diagonalization of a complex valued 
matrix which becomes computationally prohibitive inside a numerical integration. Fortunately, superconductivity is 
effected by phonon exchange between electron states which arc close to the Fermi energy and we can approximate the 
dispersion by a power series expansion about kp at each point on the Fermi surface and retain only the linear term. 
This is because the maximum phonon energy is much less than the characteristic electronic energy scale and because 
there are no strong features in the density of states close to the Fermi energy. We therefore model the electronic 
dispersion in the second and third bands using a simpler relation, 

(k-kF)^\, (1) 



where is the Fermi velocity on band i as input and 

perpendicular to the Fermi surface at k. 
Using the Anderson and Gold matrix and averaging over the Fermi surface we find a value of Vp' = —1.367 x 

10^ m s^^ < and Vp'^ = 1.055 x 10^ m s~^. The second band is hole-like and the third band electron-like. The 
location of the Fermi surface in fc-space is somewhat complicated but we can make some simplifying assumptions. 
The second band is approximated well with a sphere and we use a constant kp for Q . The third band is more 
complicated as shown in Fig. 3. The Fermi surface consists of "pipes" around the edges of the Brillouin zone. Details 
of the Fermi surface parameterization used in the calculations are given in the Appendix, section VII. 



k — kp 



is the magnitude of the component of k — kp 



III. ELECTRON-PHONON MODEL 



Use of the experimental data ensures that the important consequences of the electron-electron interactions at the 
low frequencies of interest here are included in the location of the Fermi surface and in the Fermi velocities for the 
second and third band carriers. As a result the Hamiltonian contains a pure band structure part, Ti-bandi ^-nd the 
electron-phonon Hamiltonian, Tiei-ph- The Hamiltonian is given by 

n = Hband + Hel-ph= ^ k u ^„^kua + Q ^ + 4>^yl+g+Q u a^k.' a' (2) 

k z^— 2,3 cr k ly u' (T qQX 

In Tiband, the electron energies, are given by our parameterization of the Anderson and Gold model eq. (1), 

a denotes the electron spin, ly denotes the electronic band. In "Hei-ph the and at-^^ are the usual creation and 
destruction operators for phonons and A is the branch index for the phonon spectrum. The electron-phonon matrix 
element is 



9iu.:MQ.X) = -^^^^ZI^ZI . (,-+q) T/(,,,)(g + g), (3) 

where the phonon spectrum, hw^\, and the phonon wavefunctions, f}^\, arc given by a force constant model for the 
phonon dispersion based on work done by Cowley V(yy')(g) is the Fourier transform of the screened electron- 
ion pseudopotential describing scattering between the two electron bands in the problem. We modeled the three 
different electron-phonon matrix elements, ff(22)(9i Qj 5(23) (9i 5(32) ("f? Qi ^) a-nd 5(33) ((f, Q, A), by a common 

pseudopotential function"^^ and scaling factors, so that 5(i/i/') (<f, Q, A) = ai^^^i-^ g{q,Q,X). The values of the a(j,i^/) are 
determined by fitting the solution of the gap equation to the experimentally determined values for the two bands^"^, 
1.29 mcV and 1.38 meV, and are found to be 0(22) — 0.61 , 0(23) ~ 0.50 , and 0(33) = 0.56^^. 



IV. OPTICAL CONDUCTIVITY 



The optical conductivity is given by 



a(a;)^-^{-\xico)+^J^, (4) 




FIG. 2: (Color online) Second band in Pb calculated using the matrix by Anderson and Gold 
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FIG. 3: (Color online) Third band in Pb calculated using the matrix by Anderson and Gold'^*. The more complicated Fermi 
surface structure is illustrated by the "pipe" like tubes along the zone edges. 
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where xi'^) = (0; j)) = ~* /o°° (bX^); j(0)]) s*"* the retarded current-current correlation function. In the normal 
state the d.c. conductivity, (j{Q), is finite so that the poles in the two terms exactly cancel. Consequently x(0) = 
(jjpi/ (47re^) = ne^/m*p(, where n is the density of carriers and m*pj is the optical effective mass. Gotze and Wolfle 
reformulated this result in terms of a memory function, M{uj)^ 



where 



Mi.) = -^^^y (6) 

X(0) - x(w) 

In this formulation the cancellation of the poles at w = is accomplished by introducing the memory function, M(u;), 
which remains to be determined. The analytic properties of the correlation functions require that the real part of 
M{uj), M'{uj), is odd in lo and that the imaginary part, M"(a;), is even. There are several potential contributions to 
M{ui) from different scattering mechanisms depending on the system being considered. Gotze and Wolfle considered 
a number of mechanisms including impurity scattering and electron-phonon scattering, among others'^. 



A. Electron-Phonon Contribution 



Here we are interested in the contribution from electron-phonon scattering which is characterized by the strength 
of the coupling, g. The functional form of M{lo) is found from xi'^) to second order in g by assuming that M{lo) is 
a power series in the strength of the electron-phonon interaction, g, so that, expanding the right hand side of eq. (6) 
one finds'^ 



M,^. ^ [(([j,g];[j,g]))^-(([j,g];[.7,g])L^o] ^ ^.4) 



where 



and 



dk' 



(7) 



(8) 



^ dh ''''' 



(9) 



In the superconducting state the response of the supcrcurrent to the external field leads to a delta-function in 
17(0;) at (jj = in the real part and a pole in the imaginary part. The weight of the pole is given by c^/ (47rA^) (see 
Tinkham'^^), where A is the penetration depth. Using experimental values for Pb, x(0) in the superconducting state 
is reduced from the normal state value by ~ 0.026. This leads to a negligible quantitative effect for the frequencies of 
interest here, lo > 2A, and therefore we will ignore this. 

In the following we drop the band index to simplify the notation. The correlation function in the superconducting 
state to lowest order in the strength of the electron-phonon interaction, 0{g^), is 



(10) 



where 



k.q.X 
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and 



- 4)(1 + n,^) - 4(1 - 



1 k+q 



k.q.X 



(1- W(l-/.-: 



■fkf k+q 



The and v^: are the usual BCS coherence factors. M{k,q^\) is the electron-phonon mtcraction defined earlier 
together with the derivatives of the quasiparticle spectrum in the current operator. 



M{k,q,X) 



Noh^ 
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k-k' 
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dip 
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(11) 



({[j, H] ; [j, H]))[^^ has the character of a particle-hole contribution and smoothly evolves into the normal state result 
derived by Allen and by Gotze and Wolfle"^. It is exponentially small for temperatures much smaller than the gap and 
is negligible at the temperature at which Farnworth and Timusk took their data on lead in the superconducting phase, 
0.35 K. On the other hand the combination of coherence factors and the particle-particle or hole-hole distribution 
functions ensures that {{[j , H] ; [j , H]))^^"^ vanishes as the gap vanishes. This term gives the conductivity at finite 
frequencies at low temperatures much less than Tc- Taking the limit as T ^ 0, the real and imaginary parts of Ms{i-o) 
become 



M^"{lo) 



[ k k+q k k+q k k k+q k+q 



k.q 
X 



\M{k,q,X)\^ 
X(0) 
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x(o) 



k+q 



^'T+EUq + E^ 



(12) 



B. Absorptivity 

The absorptivity is proportional to the real part of the surface impedance, (1 + (47ri/ci;) cr([j)) . As a minimal 
model we ignore the difference between the values of the gaps on the Fermi surfaces associated with the second and 
third bands and take Ag = 1.37 meV. This is an average of the two values found by Farnworth and Timusk We 
will therefore get three contributions to the memory function, M{u}), and to the absorptivity. These three are from 
intra-band scattering within the second band, M(22)(w), intra-band scattering within the third band, M(33)(a'), and 
inter-band scattering between the second and the third band, M(-23)(w). Since the frequencies of interest are much 
smaller than the plasma frequency, iVpi, the expression for the absorptivity A{uj) can be greatly simplified. 
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(13) 
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FIG. 4: (Color Online) M'^22){^) ^- 

The full line is M^^^-^ (a;) in the superconducting state while the dashed line is M^^^^{uS] 

in the normal state. 



Here we have also assumed that 1 + ^"^J"-* > 1 and < 1 so that square root can be expanded. ^^J"-* evaluated 

at = is the enhancement of the optical effective mass which is similar to the enhancement of the effective mass 
due to the electron-phonon interaction in the specific heat. 

The -M(22)j -^(23) l'^)' ^^^^ -^-^(33) contributions in the normal and superconducting states to order are compared 
in Figs. 4, 5, and 6. The different magnitudes of the memory functions Af(y)((jj). reflects their dependence on the 
surface area of the different parts of the Fermi surface involved in the scattering processes and on the size of the phase 
space for phonons which connect different parts of the Fermi surface. This dominates any effects the density of states 
on the different Fermi surfaces might have due to the cancellation of the Fermi velocity factors in the density of states 
with those from the current operator matrix elements. 

The dominant contribution comes from -^^'33-) (y>) which involves the scattering between the 36 "pipes" which make 
up the third band. Although the geometry of the third band suggests that quasi-one dimensional "nesting" , scattering 
between parallel pipes with a single phonon momentum, could be important, it is found not to dominate -^^'33^(0;). 

There are = 648 different combinations of pipes and of those, pipes whose mid-points can be connected by 
momentum vectors Q\ = 27r/a (1, 1/2, 1/2) and Q-j = 27r/a (3/2, 0, 0) lead to contributions much larger than other 
combinations by about a factor of ^ 10. There are only approximately 100 of these combinations and their contribution 
does not dominate M^33^(a;). 

In contrast to the intra-band scattering events on the second band, there are many more phonons that connect 
the more complicated surface of the third band. This and the larger surface area of the third zone Fermi surface is 
responsible for the big difference in magnitude of the two contributions. 

The real parts of the are related to the imaginary parts through a Kramers-Kronig relation, M'{lo) = 

^ dz , so that M'{lo) can, in principle, be calculated directly from M" (lo). Consequently the magnitudes of 
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FIG. 5: (Color Online) M('23) (a;) vs. 
in the normal state. 



The full line is Mfag) (uj) in the superconducting state while the dashed line is Mf^'^Juj) 



the real parts of the Afj^ij) have a similar relation to one another as the imaginary parts, hm^j^o M'{lo)/u) corresponds 
to the optical effective mass enhancement due to the electron-phonon interaction, similar to the enhancement seen in 
the specific heat. M'(uj)/uj in the superconducting and normal states are smooth functions of w which initially increase 
from their = values and reach maximum values at approximately 3A in the normal state and approximately 5A 
in the superconducting state before falling off monotonically with increasing lo. Beyond lo = 6 A they are almost equal 
in magnitude. 



C. Strong-Coupling Features 



In order to analyze the results further, we follow Farnworth and Timusk and take the derivative of A[lij) 
respect to lo. This reveals the strong-coupling electron-phonon features. 



with 



dAsjuj) 
doj 



dAN{uj) dM 
oc 



S/' 



dio 



dio 



doj 



(14) 



The strong-coupling features come almost exclusively from the first term on the right-hand side of the equation. In 
figures 7, 8, 9, and 10 we compare the contributions from the different scattering processes to dM^"{uj)/dijj. In each 
case the peaks at ~ 35 cm~^ and ~ 70 cm~^, present in the data of Farnworth and Timusk, are recovered. These 
peaks come from features in the phonon density of states and are associated with a wide range of wavenumbers. The 
M^^"(cj) in the normal state is featureless above 4A and subtracting from the superconducting M^-'J-^{uj) does not 
introduce features not already present. 
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FIG. 6: (Color Online) M'f'^^^{Lu) vs. uj. The dashed line is Mjlg-, (a;) in the superconducting state while the full line is M^^'^{lu) 
in the normal state. 

M^33j(a') is the dominant contribution to M^"{lli) and its derivative strongly resembles the data in Farnworth and 

Timusk^'*. The two prominent peaks at C and G closely resemble the peaks at the same energies in Fig. 9 from 
the point of view of both height and width. The features at B, D and E are also seen in Fig. 9 and the relative 
magnitude of the calculated dM" /dio at E and at C and G is the same as in the data. The intcrband scattering 
contribution, dM'fi^^^/duj, enhances equally the magnitudes of the peaks at C and G but shows none of the other 
features. These differences arise from the shape of the two contributions to the Fermi surface from the two bands. 
Scattering between different parts of the third band Fermi surface introduces all the phonon modes and as a result 
it is the dominant contribution and closely resembles the phonon density of states. By contrast scattering within 
the second band constrains the phase space for phonons which conserve energy and momentum. The magnitude of 
the higher energy peak in dM'^22)i^)/^^ much smaller than that of the low energy peak, whereas the peaks are 
of equal magnitude in (iA/j23^(w)/da; and dM'^^^-^{uj) / du. The high energy peak comes from phonons whose energy is 

8.48 meV. In the case of M22(ij-') the magnitude of the wavenumbers of the phonons with this energy lie between 
1.33 A^^ and 1.43 A^^. By contrast the corresponding phonon wavenumbers in the case of M33(a;) are grouped in 
sets with \q\ - 0.78 A"\ ~ 0.8 A"\ 0.86 A"\ and ~ 0.88 A"^ In Fig. 10 we show the different contributions to 
dW'/dijj. 

In Fig. 11 we plot the calculated difference of the derivatives of As{u}) and Ajv(a;) frequency shifted downwards 
by 2A for direct comparison with Farnworth and Timusk's data. We have also shifted Farnworth and Timusk's data 
by 0.9 cm~^ which takes account of the slightly different values of A used by them and us. Many of the features of 
the experimental data are reproduced: the shoulder B (~ 27 cm~^) is enhanced in our calculation, the peaks at C 
(~ 36 cm~^) and at G (~ 67 cm~^) are recovered and are equal in magnitude as in the data, and the small feature 
at E (~ 52 cm"-'^) is reproduced. The shoulder at D (~ 50 cm~^) is much weaker in the calculation compared to the 
data and the dip at ~ 56 cm~^ between E and F is missing. These discrepancies in the dM^"/duj are probably due 
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FIG. 7: (Color Online) The difference between d" M (^22) {'^) / duo (full line) in the superconducting and normal states. 



to our parameterization of the Fermi surface. As we pointed out earlier, the feature at A arises from contributions 
which are fourth order in the electron-phonon coupling and so are naturally absent from our calculation. 

The contribution to cr{Lo) from the Holstein process in Pb was investigated many years ago by Allen^^. He used the 
isotropic approximation familiar from the conventional Eliashbcrg formalism in which momentum conservation is not 
enforced in the clcctron-phonon scattering. In this approach the details of the bandstructurc arc completely absent. 
He introduced a transport version of the a^F(uj), al^F(uj), in which the high frequency components were reduced to 
allow for the suppression of forward scattering in the matrix element. He used this to evaluated a frequency dependent 
relaxation time. This amounts to putting in by hand the strong-coupling features which we have calculated by taking 
the derivatives of the Afj^ (cL;)'s. Allen used the isotropic approximation to calculate this contribution to (t{lo). His 

result, a1^F{ui) (Fig. 4 of AUen^^), is very similar to dM22{oj) / dui in Fig. 7. The main features at 4-5 meV and at 
8 meV are recovered. However, as he remarked, there is a discrepancy between the calculation of the absorptivity 
and the data of Joyce and Richards'^* and he offered a number of possible sources for this. Our calculation indicates 
that this discrepancy arises from the use of the isotropic approximation for the Fermi surface. 

Allen suggested that one source of the discrepancy between his calculation and the data could be the absence 
of strong coupling corrections described by a frequency dependent gap function as calculated from the Eliashbcrg 
equations. Our calculation, using a momentum and frequency independent gap, reproduces much of the fine structure 
seen in the strong-coupling features and it seems likely that the differences with the data may be due to simplifications 
in the parameterization of the Fermi surface or phonon spectrum. However the absence of structure in the gap function 
could lead to corrections. Allen referred to a formalism for including a gap with structure developed by Nam'^'^ which 
incorporates details of the scattering mechanism by using the single-particle self-energy determined by the Eliashbcrg 
equations. Apart from the lack of momentum conservation in this conventional application of these equations, which 
could be fixed by adding momentum dcpendence^'^, the single-particle lifetime is not the same as the transport 
relaxation time which is given by M"{lu) in our calculation. A more complete theory remains to be developed which 
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FIG. 8: (Color Online) The difference between d" M^23){^) /d^ (Ml line) in the superconducting and normal states. 

incorporates both a realistic Fermi surface and corrections to the constant gap approximation. 

Finally, we comment on the identification of M{uj), an "optical single-particle self-energy", with the single-particle 
self-energy, as extracted from ARPES by Hwang et al.^ in their analysis of data on Bi2Sr2CaCu208+5. In the 
electron-phonon model discussed here the M^"^-^(w)'s are quite different in form from the well-known imaginary part 
of the single-particle self-energy due to the electron-phonon interaction, T,"{E), as shown in Figure 6.16 of Mahan'^^. 
Although S(£') and the M{u;) have the same analytic requirements, the real parts are odd in frequency while imaginary 
parts are even, so that they have similar behavior at low frequencies, their finite frequency behavior is quite different. 
The imaginary part of the self-energy is constant beyond ~ 10 meV or ^ 7.3Ao whereas -^("j) (<^)'s grow monotonically 
as the phase space for quasiparticle states increases with uj. Further, in the superconducting state, a feature associated 
with a peak in the phonon density of states, which appears at ujph (say), appears at A + ujph in Y."{E) but at 2A + u}ph 
in M^"{uj). If further analysis of data on Bi2Sr2CaCu208+5 continues to support the proposed relation between M{uj) 
and S(£'), this will provide a strong constraint on the imderlying Hamiltonian. 

V. CONCLUSIONS 

Using the memory function formalism we have demonstrated the importance of including a realistic Fermi surface 
in Pb to obtain quantitative agreement with the strong-coupling features seen in the optical data. The largest 
contribution to the strong-coupling features comes from intra-third band scattering which also provides most of the 
fine structure. The phonons involved in scattering within this band are not as restricted by momentum conservation as 
those involved in intra-second band scattering because the third band contribution to the Fermi surface is spread over 
a substantial fraction of the region near the zone boundary. This seems to be the reason why the a^Fioj) determined 
from tunneling and that determined from a{uj) are very similar, whereas one might have assumed that suppression of 
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FIG. 9: (Color Online) The difference between d"M(33) (w)/da' (full line) in the superconducting and normal states. Note that 
the peak at ~ 8.5 meV has the same height as that at 5.5 meV. This reproduces an important feature of the Farnworth and 
Timusk data^''. 



forward scattering in the current matrix elements would have lead to a difference between them. Allen's investigation 
of the relation between a'^F{uj) and aj^F(w) assumed a spherical Fermi surface and his model for a^rF{uj) is very 
similar to the strong coupling features found from scattering between states on a spherical Fermi surface, dM'^^^^-^/ (Lj 
above. 

Strong-coupling features in tunneling and transport data are a useful check on microscopic models of the mech- 
anism for superconductivity. There are a number of systems in which mechanisms other than the electron-phonon 
interaction are thought to play a role. Among these are the cuprates where the contribution of the electron-phonon 
interaction remains to be clarified and in which the magnetic properties have been proposed as the origin of the 
d-wave superconductivity. Analysis of data so far finds no evidence of electron-phonon strong coupling features'^^, 
although Maksimov et al.'^^ have proposed a model for the electron-phonon interaction with which they have analyzed 
ARPES data on Bi2Sr2CaCu208. The data of Hwang et al.^ show that M"{uj) for Bi2Sr2CaCu208+5 has the same 
monotonic frequency dependence as was found in current work on Pb at low temperatures so that its derivative may 
reveal strong-coupling features associated with a magnetic mechanism. Analysis of data on the cuprates to date has 
employed the same formalism as in the electron-phonon problem with the atrF{uj) function replaced with a spectral 
density for spin fluctuations. The weak coupling approximation introduced by Gotze and Wolflc"^ would have to 
be extended to calculate the contributions to M(a;) from electron-phonon or impurity scattering enhanced by spin 
fluctuations or from Umklapp scattering due to interactions among the charge carriers. Riseborough^^ has previously 
calculated the contribution from Umklapp scattering due to spin fluctuations in the random phase approximation 
using the Kubo formalism. A memory function approach extended to deal with intermediate to strong coupling could 
provide a method to test microscopic models for spin fluctuation mechanisms for superconductivity in the cuprates 
while taking account of bandstructure and momentum conservation. 
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FIG. 10: (Color Online) Contributions to dAs{ij)/dLo — dAN{Lo) /dcj plotted against lo — 2A in units of cm ^(fuU line) from 
the three types of scattering processes. 
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FIG. 12: The Fermi momentum as a function of angle /3 for the two types of third band pipes. 

VII. APPENDIX: FERMI SURFACE PARAMETERIZATION 

Although diagonalizing the Anderson and Gold 8x8 matrix is computationally not very challenging in itself, we will 
need to calculate its eigenvalues many times since we will integrate over the electron momentum. It turns out that 
using a diagonalization routine directly inside the integration is computationally not feasible. We therefore had to 
simplify the model. 

As pointed out earlier, we will require two of the four bands only. The second band Fermi surface is close to 
spherical which led us to set the Fermi momentum in that band to a constant, Vp'' = —1.367 x 10^ m s~^. The second 
band is hole-like. 

The third band Fermi surface consists of "pipes" which we approximated by dented cylinders. There are two types 
of pipes and we chose to parametrize these separately. In Fig. 12 we show the angular dependence of the Fermi 
momentum on each of the two pipes. The dispersion was linearized and the Fermi velocity we calculated for both 
types is = 1.055 x 10^ m s~^. The third band is electron-like. 



